Statistical properties of 12.2 GHz methanol masers associated 
with a complete sample of 6.7 GHz methanol masers 



S. L. Breen, 1 ' 2 S. P. Ellingsen, 2 J. L. Caswell, 1 J. A. Green, 1 G. A. Fuller, 3 M. A. 

Voronkov, 1 L. J. Quinn, 3 A. Avison 3 
1 CSIRO Astronomy and Space Science, PO Box 76, Epping, NSW 1710, Australia 
2 School of Mathematics and Physics, University of Tasmania, Private Bag 37, Hobart, 

Tasmania 7001, Australia 
3 Jodrell Bank Centre for Astrophysics, Alan Turing Building, School of Physics and 
Astronomy, University of Manchester, Manchester M13 9PL 

Shari . BreenOcsiro . au 



ABSTRACT 

We present definitive detection statistics for 12.2 GHz methanol masers towards a complete 
sample of 6.7 GHz methanol masers detected in the Methanol Multibeam survey south of decli- 
nation -20°. In total, we detect 250 12.2 GHz methanol masers towards 580 6.7 GHz methanol 
masers. This equates to a detection rate of 43.1 %, which is lower than that of previous signif- 
icant searches of comparable sensitivity. Both the velocity ranges and the flux densities of the 
target 6.7 GHz sources surpass that of their 12.2 GHz companion in almost all cases. 80 % of the 
detected 12.2 GHz methanol maser peaks are coincident in velocity with the 6.7 GHz maser peak. 
Our data support an evolutionary scenario whereby the 12.2 GHz sources are associated with a 
somewhat later evolutionary stage than the 6.7 GHz sources devoid of this transition. Further- 
more, we find that the 6.7 GHz and 12.2 GHz methanol sources increase in luminosity as they 
evolve. In addition to this, evidence for an increase in velocity range with evolution is presented. 
This implies that it is not only the luminosity, but also the volume of gas conducive to the different 
maser transitions, that increases as the sources evolve. Comparison with GLIMPSE mid-infrared 
sources has revealed a coincidence rate between the locations of the 6.7 GHz methanol masers 
and GLIMPSE point sources similar to that achieved in previous studies. Overall, the proper- 
ties of the GLIMPSE sources with and without 12.2 GHz counterparts are similar. There is a 
higher 12.2 GHz detection rate towards those 6.7 GHz methanol masers that are coincident with 
extended green objects. 

Subject headings: masers — stars: formation — ISM: molecules — radio lines: ISM 



1. Introduction 

Masers of water, methanol and hydroxyl are 
considered excellent probes of high-mass star for- 
mation regions. They are common, bright and 
are observed during the phase where the forming 
stars are still deeply embedded in their natal ma- 
terial. The different maser species require slightly 
different physical co nditions for their produc- 
tion and survival (e.g. Cragg. Sobolev k, Godfrey! 
2002), which have led to proposals which base 



high-mass star formation evolutionary schemes 
on the presence or absence of the different maser 
transitions (e.g. Ellingsen et al.l 2007t Breen et al 



2010at iFontani. Cesaroni fe Furuvall2010l ) 



Previous searches for 12.2 GHz methanol 
mase rs have typi c ally t argeted UCHn reg ions 



(e.g. iBatrla et al.l 119871 : iNorris et al.l 119871 ) or 



6.7 GHz methanol masers that were detected to- 
wards UCHn r egions, OH masers or IRA S selected 
regions (e.g. iGavlard. MacLeod fc van der Walt 
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1994t ICaswell et all Il995bt iBlaszkiewicz fe Kusl 
2004; Br een et al.ll2010a[ ). OH masers have been 
noted on many occasions as being associated 
with a generally later stage of star formation 
than other m aser sp ecies (e.g. iForster fe Caswell 
19891: ICaswelJfl99"7t) . Furthermore, those IRAS 



promising resul t s (e. 



Chambers et al. 200 



Cvganowski et al 



Rathbornc et al 



2008; 
2010h . 



selected searches were targeted towards sources 
that h ad colours similar to UCHi i regions based 
on the IWood fe Churchwelll ( 19891 ) criteria. Each 
of the aforementioned 12.2 GHz methanol maser 
searches have therefore been biased towards 
sources of a generally later stage of high-mass 
star formation since we know that many 6.7 GHz 
methanol masers are devoid of one o r both of as- 
sociated OH masers dCaswelll fl996l) and UCHn 
regions (jWalsh et al.lll998l ). 

Breen et al. (2010a) used their 12.2 GHz methanol 



maser observations (together with complemen- 
tary data) to investigate the relative evolutionary 
stages of the high-mass star formation that the 
6.7 and 12.2 GHz methanol maser transitions were 



tracing. From this analysis I Breen et al.l (|2010af ) 
concluded that 12.2 GHz methanol masers were 
present towards the second half of the 6.7 GHz 
methanol maser lifetime. They additionally found 
that both methanol transitions appeared to in- 
crease in luminosity as the sources evolved, al- 
though at a slower rate and by a lesser amount 
in the case of the 12.2 GHz sources. Although 
comprising a relatively small and biased sample of 
sources, these results are very encouraging in the 
overall goal to use the common maser species as 
evolutionary tracers of high-mass star formation. 

The idea that the presence, absence or inten- 
sity of certain properties associated with high- 
mass star formation regions can give an indica- 
tion on the relative evolutionary stage of star for- 
mation is not a new one. In particular, from 
soon after the discovery of interstellar masers, 
it has been suggested that these masers are as- 
sociated with a dist inct evolutionary phase of 
star formation (e g. iLo. Burke fc Haschick 19751: 



Genzel fc Downeslll977l) and since this time, many 



maser studies have found evidence that this is the 



case (e.g. Forster & Caswelll Il989: 


Caswelll 19971 


Walsh et al. Il998l Ellingsenl 12006 


: iBreen et al. 


2010al Fontani. Cesaroni & Furuva 


20101 Evolu- 



tionary scenarios based on combinations of specific 
infrared, submillimetre and millimetre radio con- 
tinuum properties are also popular and show very 



Thermal lines have also proven to be excellent 
evolutionary probes, especially since, in addition 
to the presence or absence of specific transitions, 
line properties can be used to further pinpoint 
the source's evolution (e.g . IPurcell et al ] 120061 
20091 : iLongmore et al.ll2007h . The implied evolu- 
tionary stage that is deduced by these different 
methods appear to be in e x cellen t agreement. For 
example, ILongmore et al. (12007ft found that the 
linewidths of the detected ammonia emission in- 
creased with evolution. Comparison between am- 
monia linewidths and m ethanol maser luminosi- 
ties by fWu et al.l (l2010h show that the lower lu- 
minosity 6.7 GHz methanol masers are associated 
with smaller ammonia linewidths and therefore 
correspond to the y ounger sources. Thi s corrobo- 
rates the findings of Breen et al. ( 2010al ). that the 
6.7 GHz methanol masers increase in luminosity 
as they evolve. 

In this paper we present the statistical prop- 
erties of the 12.2 GHz methanol maser emission 
detected towards a complete sample of 6.7 GHz 
methanol masers. In particular, we t est the con- 
clusions drawn in IBreen et al.l (|2010a|) on a much 
larger and less biased sample of sources. Through- 
out this work we place a particular emphasis 
on investigating the properties of the different 
maser species at different evolutionary stages. The 
6.7 GHz methanol maser targets comprise a com- 
plete sample of almost 600 methanol masers de- 
tected in the Methanol Multibeam Survey (de- 
scribed in Section [P]) . In a subsequent publica- 
tion we will be presenting the catalog of detected 
12.2 GHz methanol masers, together with discus- 
sion of individual source properties that can be 
best addressed by looking at their specific data. 
This publication order is necessary since many of 
our target sources are yet to be published. 

1.1. The Methanol Multibeam (MMB) 
Survey 

The Methanol Multibeam (MMB) Survey for 
6.7 GHz methanol masers is the largest complete 
survey for 6.7 GHz methanol masers ever under- 
taken. The complete survey aims to cover the full 
Galactic Plane within b = ±2°. A thorough de- 
scription of the survey parameters and techniques 
can be found in Green et al. I (120091 ). 
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The southern hemisphere component of the 
MMB survey has detected nearly 1000 6.7 GHz 
methanol masers (to a detection limit of ~0.7 Jy), 
the majority of which have now been followed 
up with interferometric observations in order to 
obtain accurate source positions. Survey results 
from the Galactic longitude range 345° to 20° 
have already bee n released ( Caswell et al. 2010t 
Green et al. 2010l ) and results from the range 186 
to 345° are forthcoming. 

2. Observations and data reduction 
2.1. 6.7 GHz MMB maser targets 

12.2 GHz follow up observations towards 593 
6.7 GHz methanol masers detected in the MMB 
survey were carried out with the Parkes radio 
telescope during 2008 June and December. The 
sources searched correspond to all of the MMB de- 
tections that lie south of declination -20 degrees. 
This declination cutoff was imposed for several 
reasons, foremost of which was that at the time 
of the 12.2 GHz observations, accurate positions 
were chiefly limited to those sources that fell in 
this range. Since the positional uncertainties of 
sources detected in the initial survey were com- 
parable to the half power beam width (HPBW) 
of the Parkes radio telescope at 12.2 GHz, follow 
up observations towards the sources that lie north 
of declination -20 degrees have been delayed until 
accurate positions are known. 

As part of the MMB observing strategy, once 
accurate positions for the survey detections were 
obtained (chiefly using the Australia Telescope 
Compact Array, achieving positional accuracies 
of ^0.4 arcsec), pointed observati ons of each 



sourc e were completed with Parkes (jGreen et al 
2009). The resulting spectra have typical rms 
noise levels of ^0.07 Jy and spectral resolutions 
of 0.11 km s _1 . It is with these spectra that we 
compare our follow-up 12.2 GHz observations. 

The 6.7 GHz methanol maser properties (ve- 
locity ranges and peak, peak flux density and in- 
tegrated flux density) of the sources targeted for 
12.2 GHz emission have been determined for the 
present analysis independently of those listed in 
the surv ey publications that have already bee n 
released ( Caswell et al. 2010t Green et al. 2010l) . 
However large discrepancies between values are 
not expected. For some sources in regions of confu- 



sion, the interferometric data had to be inspected 
to determine which spectral features were associ- 
ated with which position. 

2.2. 12.2 GHz observations 

The search for 12.2 GHz methanol masers were 
made with the Ku-band receiver on the Parkes 
64 m radio telesc ope, using an identic al setup to 



that described in iBreen et al.l ( 2010al ). This re- 



ceiver detected two orthogonal linear polarizations 
and had typical system equivalent flux densities 
of 205 and 225 Jy for the respective polarizations 
throughout the observations in 2008 June, and 
were slightly higher at 220 and 240 Jy in 2008 
December. The Parkes multibeam correlator was 
configured to record 8192 channels over 16 MHz 
for each of the recorded linearly polarized signals. 
This configuration yielded a usable velocity cov- 
erage of ^290 km s _1 and a spectral resolution of 
0.08 km s , after Hanning smoothing of the data. 
The Parkes radio telescope has rms pointing errors 
of ~10 arcsec and at 12.2 GHz the telescope has 
a HPBW of 1.9 arcmin. For two sources, no us- 
able data was obtained in either the 2008 June or 
December observations. Instead, data from sub- 
sequent observations made during 2010 March are 
used. 

The observation strategy was crafted such that 
the majority of the target 6.7 GHz methanol 
masers had a total of 10 minutes of integration 
time at 12.2 GHz over the two main observing 
epochs. This integration time yields typical rms 
noise limits of ~0.11 Jy (corresponding to a 5- 
cr detection limit of 0.55 Jy). A mixture of 10 
and 5 minute integrations were taken at the first 
epoch. Sources observed for only 5 minutes (dur- 
ing the first epoch) which did not yield a detection 
were typically observed for an additional 5 minutes 
during the second epoch. Confident detections 
were not routinely repeated, regardless of their 
integration times. For 55 of the target 6.7 GHz 
methanol masers, the 12.2 GHz methanol maser 
d ata used in th i s analy sis is the same as published 
in lBreen etall l|2010af ). 

All sources were observed at a fixed frequency of 
12178 MHz (i.e. with no Doppler tracking), which 
alleviated the requirement for a unique reference 
observation to be made for each of the source po- 
sitions. Instead, all the spectra collected each day 
were combined to form a sensitive reference band- 
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pass by taking the median value of each channel. 
Since the target sources were spread over a large 
Galactic longitude range, sources were detected at 
different frequencies within the bandpass and this 
combined with the narrow nature of the sources 
means that this technique works very well, even 
in the case where a high number of strong detec- 
tions are expected. The noise contribution to the 
quotient spectrum from such a reference is negli- 
gible. 

Data were reduced using ASAP (ATNF Spec- 
tral Analysis Package). Alignment of veloc- 
ity channels was carried out during processing. 
The adopted rest frequency wa s 12.178597 GHz 
(jMuller. Menten fc Madeil |2004| ). Absolute flux 
density calibration was achieved by observing PKS 
B1934-638 each day which has an assumed flu x 
density of 1.825 Jy at 12178-MHz (|Sauld l2003h . 
Calibration observations were very stable during 
each session. 

Sources that were observed at both epochs were 
first reduced individually and searched for de- 
tections and then the data were averaged and 
searched once more in order to pick up any sta- 
ble but weak emission. Only three sources were 
identified in this manner (out of 208 sources that 
were observed at both epochs with no detection at 
either epoch). Given the luxury of knowing the ve- 
locity of the 6.7 GHz maser emission, we have been 
able to readily and reliably identify sources that 
are much weaker than the 5-er detection limits. A 
number of weak sources can confidently been re- 
garded as genuine at 3- to 4-er levels. Marginal 
detections were considered along with their asso- 
ciated 6.7 GHz spectrum, and were assessed on an 
individual case basis. 

3. Results 

Observations conducted with the Parkes radio 
telescope over two epochs resulted in the detec- 
tion of 250 12.2 GHz methanol masers towards 580 
6.7 GHz methanol masers detected in the MMB 
survey, a detection rate of 43 %. A search of the 
literature indicates that 148 of the 250 12.2 GHz 
methanol masers that we detect are new discov- 
eries. Of the 102 12.2 GHz sources that have 
been previously discove red, the major i ty of th ese 
were first detected by I Caswell et al. ( 1995bl ) or 
Breen et"al~l ( 2010al ). but others were first detected 



bv: IGavlard. MacLeod fe van der Waltl (ll994Tl; 



iKoo et al l (Il988l): ICatarzi. Moscadelli k Panella 
(|l993h : lKemball. Gavlard fe Nicolsonl(ll988[):lBatrla et al 
(ll987li: lCaswell et all fll993h: iNorris et all (Il987h:~ 



MacLeod. Gavlard fe Kemballl (Il993h . The ob 



or 

servations targeted a total of 593 sources south of 
our declination cutoff of -20 degrees but 13 sources 
have been removed from the detection statistics, 
chiefly due to our inability to distinguish them 
from nearby sources (a consequence of our single 
dish observations). Ten of these 13 sources are 
located close to the Galactic centre. 

4. Discussion 



4.1. Detection statistics 

Our detection rate of 43 % is significantly lower 
than previous searches of comparable sensitiv- 
ity. For example, the detec tion rates found in the 
larg e, sensitive surveys of ICaswell et~aT1 (|l995bh 
and iBreen et al. ( 2010al ) were 55 % and 60 % re- 
spectively. In order to determine the reasons for 
this apparent detection rate discrepancy, we must 
first examine the natu re of the target sources used 



in previous searches. I Caswell et al. I (Il995bh 



ducted their search towards 6.7 GHz methanol 
masers that had been detected towar ds OH maser 
targets , while the observations of IBreen et al 



(|2010al ) targeted 6.7 GHz methanol masers that 
had been detected t oward s I RAS selected regions 
(IWalsh et al.l Il997l Il998l ). iBlaszkiewicz fe Kusl 
(J2004J) searched 261 6.7 GHz methanol masers 
and achieved a detection rate of 19 %. Their 
survey was of much lower sensitivity (account- 
ing for the lower detection rate) and was also 
primarily directed towards 6.7 GHz methanol 
masers detected towards IRAS se l ected targets. 
IGavlard. MacLeod fe van der Walt ( 1994 ) carried 
out a search for 12.2 GHz methanol masers to- 
wards infrared selected 6.7 GHz targets with a 
similarly low sensitivity, and achieved a detection 
rate of 27 %. 

Given the high sensitivity and large scale of 
our observations, we can only directly c ompare 
our results with those of Caswell et al. ( 1995bl) 
and IBreen et al.l (|2010al ); however, we are able to 
make some comparison with less sensitive searches 
by taking detection limits into account. For ex- 
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ample, iGavlard. MacLeod fc van der Waltl (|1994l ) 
had a sensitivity limit of ~10 Jy, if our detection 
limit was at this level our detection rate would be 
slightly lower than 10 % (compared to their detec- 
tion rate of 27 %). If our dete ction limit had been 
4 Jy, l ike that of the search bv lBlaszkiewicz fc Kusl 
(|2004h . our detection rate would have dropped to 
16 % (compared to their detection rate of 19 %). It 
is therefore evident that even by taking detection 
limits of other surveys into account, our detection 
rate is lower than that of previous searches. 

T he large sea r ches o f ICaswell et all (|l995bh 
and Breen et al. ( 2010al ) were targeted towards 
6.7 GHz methanol masers that were detected in 
observa tions of comparabl e sensit ivity (|Caswell et al 
(jl995bf ) and IWalsh et all (|l997l) respectively) to 
that of our target 6.7 GHz sources. This means 
that the difference in 12.2 GHz detection statis- 
tics cannot be accounted for by the sensitivity in 
the searches of the target 6.7 GHz sources. The 
marked drop in detection rate from our search 
compared to similar searches can therefore most 
easily be accounted for by the absence of 6.7 GHz 
methanol maser target biases that were present 
in many previous searches. Our 6.7 GHz MMB 
targets were detected in a sensitive and unbiased 
manner. In contrast, all of the previous searches 
were subject to different levels of evolutionary bi- 
ases, since both IRAS infrared colour and OH 
targeted searches would miss the majority of very 
young methanol masers. Add to this the fact that 
12.2 GHz methanol masers are expected to be 
present at a time somewh at later than the on set of 
the 6.7 GHz emission (e.g. lBreen et al.l2010al) . and 
it is not surprising that our detection rate would 
be low er than such searches. iBlaszkiewicz fc Kusl 
(|2004l ) suggested that 12.2 GHz detection rates de- 
pended mainly on sensitivity; however, our search 
shows that the role of biases is perhaps an equally 
important factor. 

4.2. Galactic distribution of methanol 



Pestalozzi et al. ( 20051 ) compiled a catalog of 



the 6.7 GHz methanol masers that were known 
prior to 2004. A total of 519 6.7 GHz sources were 
listed throughout the Galaxy and were detected in 
both targeted and blind searches (not all of which 
had been searched for 12.2 GHz counterparts). 
The distribution of the 6.7 GHz sources was found 



by the authors to be non-random throughout the 
Galactic plane with a particularly high concen- 
tration betwee n 20 to 50° long i tude and 330 to 
340° longitude. IPestalozzi et all (|2005l ) also found 
that the latitude distribution of sources across 
the Galactic plane had a FWHM of 0.5° and 
commented that this could be attributed to the 
fact that blind searches had been limited to low 
Galactic latitudes. However, a comparable nar- 
row latitude distribution is seen in the MMB 
(|Caswell et al.ll2010i: [Green et alj|2010h suggesting 
it is inherent to the population rather than a prod- 
uct of observational bias. 

The Galactic distribution of the 6.7 GHz 
methanol masers detected in the MMB survey 
is beyond the scope of this study and we only 
include preliminary findings which allow us to 
comment on the distribution of the associated 
12.2 GHz sources. The sample of MMB 6.7 GHz 
methanol masers south of -20° (corresponding to 
longitude ranges of 290° through 360° to 10°) does 
not include the densely populated 20° to 50° longi- 



tude r ange that is covered by the IPestalozzi et al 



(|2005l ) compilation, but encompasses a larger 
group of sources. In the overlapping ranges, the 
distribution is similar with a significant density of 
sources in the 330 to 340° longitude range, and a 
distinct lack of sources between longitudes of 250 
to 290°. 

A comparison of the detection rate of 12.2 GHz 
methanol masers split into 5° longitude bins shows 
that the detection rate is somewhat dependent 
on Galactic longitude range (shown in Fig. [TJ. 
Sources with longitudes less than 290° or more 
than 10° were omitted from this investigation as 
these 5° bins contained fewer than 5 sources, mak- 
ing the statistics meaningless. 

In Fig. [U the the diagonally shaded bars show 
the range of detection rates in each longitude range 
that would not be significantly different from the 
mean detection rate. Therefore, in order for the 
detection rate to be statistically significantly dif- 
ferent from the mean, it must fall outside of these 
diagonally shaded bars. For the 335 to 340° range, 
the detection rate differs from the mean detec- 
tion rate significantly (p-value=0.04 from a chi- 
squared test) and shows a relative overabundance 
of 12.2 GHz sources. This may indicate that a 
burst of high-mass star formation occurred in this 
region at a time which now causes an abnormally 
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high fraction of sources to be at similar evolution- 
ary stages. 

An apparent relative under-abundance of 12.2 GHz 
detections at longitudes between 270 to 305° was 
further investigated and compared to a region of 
the same size, 315 to 350°. The detection rate in 
the first region is 25 % and in the second the de- 
tection rate almost doubles to a rate of 48 % and 
neither rate differs significantly from the mean. 
However, the detection rates in these regions are 
statistically different from one another (tested by 
calculating the expected number of detections in 
270 to 305° using the detection rate of the 315 to 
350° region and carrying out a chi-squared test). 
In the 270 to 305° region, there is also a distinctly 
lower density of 6.7 GHz methanol masers than 
in other parts of the Galaxy. In the Large Magel- 
lanic Cloud (LMC), the heavy metal and complex 
molecule abundances are a likely explanation for 
the lack of 6.7 GHz methanol masers, which are 
undcrabundant by a factor of ~4-5 compared to 
our Galaxy afte r accounting for th e different star 
formation rates ([Green et al. 2008 ). 

As evident in IPestalozzi et al. ( 20051 ) and in- 
spection of the distribution of the MMB sources, 
relatively few 6.7 GHz methanol masers reside 
in the longitude range 270 to 305°. Consider- 
ing the explanation of the underabundance of 
methanol maser sources located in the LMC, it 
is possible that the underabundance of 6.7 GHz 
methanol masers in the 270 to 305° region could 
be attributed to low metallicity, which if it is 
the cause, may also be responsible for influencing 
the 12.2 G Hz detection rate (whic h is lower than 



expected). lEllingsen et al. ( 201dl ) carried out 



search for 12.2 GHz methanol maser emission to- 
wards the four known 6.7 GHz methanol masers 
in the LMC, detecting one source. Comparison 
of the 12.2 GHz detection rate in the LMC com- 
pared to that in the Galactic plane is difficult since 
the sample size in the LMC is small and achieving 
sensitivities comparable to sources located within 

" (120101) 



Ellingsen et al 



our Galaxy is difficult, 
comment that observations of two orders of mag- 
nitude greater sensitivity would be required to 
confidently determine if the sources in the LMC 
were under-luminous at 12.2 GHz compared to the 
Galactic sources. 



4-2.1. Clustering of sources exhibiting the same 
methanol transitions 

The locations of the two methanol maser tran- 
sitions were further inspected on a longitude- 
latitude plot. This plot showed an apparent clus- 
tering of the sources that exhibited masers at both 
transitions, and likewise for the sources exhibiting 
only 6.7 GHz emission. For all sources harboring 
12.2 GHz emission, the nearest neighbor spatially 
has been examined to determine if there is a higher 
proportion of sources with 12.2 GHz emission near 
other 12.2 GHz sources than should be expected 
from a random distribution. Using this method we 
find that approximately 50 % of 12.2 GHz sources 
have another source with 12.2 GHz emission as its 
nearest neighbor (on the longitude-latitude plot) 
compared to the expected proportion (based on 
detection statistics) of 43 %. The difference be- 
tween these values is not quite statistically signif- 
icant. 

The apparent clustering was further tested by 
investigating the nearest neighbor of each of the 
12.2 GHz methanol masers on a longitude- velocity 
plot. This was done by determining the linear sep- 
aration between each of the 12.2 GHz methanol 
sources and all other sources on the plot. The dif- 
ference in longitude between source was multiplied 
by two, making the ranges on each of the axes ap- 
proximately equivalent. We find that for 51 % 
of 12.2 GHz methanol sources, the nearest source 
also exhibits 12.2 GHz emission which is statisti- 
cally significantly higher that what would be ex- 
pected from the 12.2 GHz detection rate (p-value 
of 0.05). Since we have limited our investigations 
to two dimensions in each case it is difficult to 
draw meaningful conclusions, but they do suggest 
that further, more sophisticated three dimensional 
analysis is warranted once distances to the sources 
are known. 

Clustering of sources showing 12.2 GHz methanol 
masers, combined with the non-uniform detection 
statistics along the Galactic plane (as discussed 
in Section I4.2|) . may imply that localised bursts 
of high-mass star formation have resulted in a 
number of nearby sources being at similar evo- 
lutionary stages. Another explanation may be 
that localised metallicity gradients influence the 
12.2 GHz methanol maser detection statistics. 
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4.3. Comparison of flux densities and ve- 
locity ranges 

Previous searches have shown that 12.2 GHz 
methanol maser emission is generally weaker 
than their 6.7 GHz counterparts. A large sam- 
ple of spectra demon strating this is presented by 
Caswell et al. I (Il995d) and illustrate many of the 
points made here. We, too, find that 12.2 GHz 
methanol masers are, in general, weaker than their 
6.7 GHz counterparts and that our 12.2 GHz de- 
tections are primarily at the stronger 6.7 GHz 
sites. Table [1] presents both the average and me- 
dian peak flux densities of the 6.7 GHz methanol 
masers (of the full sample as well as in the cat- 
egories of with and without 12.2 GHz counter- 
parts) and the 12.2 GHz sources that we detect. 
The average peak to peak ratio of 6.7 GHz to 
12.2 GHz methanol masers (for detected sources) 
is 18:1 and the median ratio is 4.3:1. Comparison 
of integrated flux densities shows an even greater 
difference in 6.7 and 12.2 GHz methanol maser 
strength; the average 6.7 to 12.2 GHz integrated 
flux density is 41:1 and have a median ratio of 
10:1. 

We find that the average velocity range of 
the 580 target 6.7 GHz methanol masers is 6.9 
km 8~ , and they have a median range of 5.7 
km s _ . There are only 8 instances where the 
6.7 GHz velocity range surpasses 20 km s , and 
in fact, this number only increases to 41 for veloc- 
ity ranges greater than 15 km s _1 . The detected 
12.2 GHz methanol maser sources have smaller ve- 
locity ranges; only 22 sources exceed a velocity 
range of 10 km s _1 . The average velocity range 
of the 12.2 GHz detections is 3.2 km s _1 , and 
the median range is 1.7 km s -1 . Interestingly, 



Table 1: Average and median peak flux densities of 
both the 6.7 GHz and 12.2 GHz methanol masers. 



Maser 


Peak flux density (Jy) 


Category 


Average 


Median 


All 6.7 GHz sources 


48 


5 


6.7 GHz with 12.2 GHz 


106 


15 


6.7 GHz no 12.2 GHz 


6.8 


2.5 


All 12.2 GHz sources 


23 


2.0 



there are very few instances (14 of 250) where 
a 12.2 GHz methanol counterpart is detected to- 
ward a 6.7 GHz methanol maser with a velocity 
range less than 2 km s _1 , even though 112 6.7 GHz 
methanol masers fall into this category. There- 
fore, we find a much lower 12.2 GHz detection 
rate (12.5 %) towards 6.7 GHz methanol masers 
with narrow velocity ranges (not greater than 2 
km s -1 ). Not surprisingly, these sources have an 
average peak flux density that is lower than that 
of the sources with the larger velocity ranges. 

Since the 12.2 GHz methanol masers are gen- 
erally weaker than the target 6.7 GHz methanol 
masers, one might expect that sensitivity limita- 
tions have prevented us from detecting 12.2 GHz 
emission towards the weaker 6.7 GHz methanol 
masers. Alternatively, this effect may be lim- 
ited to only the weaker features of some 6.7 GHz 
methanol masers. If these weaker 6.7 GHz features 
were located towards the minimum and/or max- 
imum velocities of the emission, then this might 
affect the recorded velocity ranges of the 12.2 GHz 
sources. However investigations presented in Sec- 
tions 233] and ElO] show that, in the majority of 
cases, sensitivity has not limited us from gaining a 
complete picture of the two methanol maser tran- 
sitions. Furthermore, given that the integrated 
flux density ratios of the two transitions are vastly 
different, this suggests that even in the cases where 
we detect both methanol transitions there is a 
larger gas volume conducive to the 6.7 GHz tran- 
sition than the 12.2 GHz transition. 

4-3.1. 6.7 and 12.2 GHz methanol maser peak ve- 
locities 

Approximately 80 % of the 12.2 GHz peak ve- 
locities correspond to the peak velocity of the as- 
sociated 6.7 GHz methanol maser. In all but one 
case, the peak velocity of the remaining 20 % of 
12.2 GHz sources corresponds to the velocity of 
another spectral feature of the 6.7 GHz source. 
In some of these cases, it is the weakest feature 
in the 6.7 GHz maser spectrum that shows emis- 
sion at 12.2 GHz, but more commonly we find the 
12.2 GHz emission at the velocity of the second 
strongest 6.7 GHz methanol maser feature. 

There are very few weak 12.2 GHz sources with 
stable peak flux densities (i.e. those that show 
variations of less than the noise level of the obser- 
vations). 236 of the 6.7 GHz methanol masers were 
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searched at both epochs and 28 of these resulted in 
a detection at only one of these, while the remain- 
ing 208 sources showed no identifiable emission at 
either epoch. When the data for these sources 
were averaged over the two epochs, only three ad- 
ditional sources were identified. This shows that 
multiple epoch searches have a much higher merit 
than one-time searches, even when the observa- 
tions are of a higher sensitivity. 

Fig- H] presents histograms that show the veloc- 
ity of the 6.7 and 12.2 GHz methanol maser peak 
with respect to the velocity range of the 6.7 GHz 
methanol maser. Sources with 6.7 GHz methanol 
maser velocity ranges of less than 2 km s _1 have 
been excluded from this plot to ensure that 
the analysis is not dominated by single-feature 
sources. Percentages of the velocity ranges are 
calculated from the most negative velocity of the 
6.7 GHz maser emission defined by the equation: 



maser is statistically significantly higher than the 
expected number (p- value 0.01). The combined 
bins between and 20 % show no significant de- 
viation from the mean. Interestingly, the two bins 
in each histogram that either achieve statistical 
significance, or close to it, all fall close to the min- 
imum velocity range. 



peak Vmin) 

percentage = 100— - — — ^> 

\ * max V rain ) 

where Y m in and V max correspond to the min- 
imum and maximum velocities of the 6.7 GHz 
methanol maser emission, and V pea k corresponds 
to the velocity of the 6.7 GHz peak and the 
12.2 GHz methanol maser peak in the top and 
bottom panels of Fig. [2] repecitvely. These his- 
tograms show that there are slightly more sources 
whose peaks lie close to either edge of the veloc- 
ity range than in the center, for both the 6.7 and 
12.2 GHz sources. 

We have compared the number of sources in 
each of the velocity range bins with the number 
of sources expected in each bin (taken to be even 
across the bins) and tested the significance with a 
chi-squared test. More 6.7 GHz methanol maser 
peaks lie close to the minimum velocity range, be- 
tween and 20 %, than the expected number (43) 
which is almost statistically significant (p-value 
of 0.07). Furthermore, there are fewer 12.2 GHz 
methanol maser peaks that lie in the to 10 % 
bin of the 6.7 GHz methanol maser velocity range. 
The number of sources that fall in this bin (14) is 
almost statistically different from the number that 
would be expected (23). However, the number of 
12.2 GHz methanol maser peaks that fall within 
10 to 20 % of the velocity range (measured from 
the minimum velocity) of the 6.7 GHz methanol 



S 




Fig. 1— Plot of the 12.2 GHz detection rates 
as a function of Galactic longitude ranges. The 
gray bars correspond to the actual detection rates 
in each longitude range and the dashed horizon- 
tal line represents the average detection rate of 
12.2 GHz sources. For the detection rate to 
be significantly different that the mean detection 
rate, they would have to fall onside the diagonally 
shaded bars (p-value=0.05). 
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Although perhaps not formally significant, a 
slight excess of sources with peak emission blue- 
shifted with respect to the systemic velocity might 
result if there are some sources where blue-shifted 
emission is enhanced as a result of an outflow, and 
the presence of a central continuum source that 
is preferentially amplified by the emission on the 
side of the source mov ing towards us (as in maser 
model s presented 
(|!997ri l. 



Sobolev. Cragg fc Godfrey 



4.4. Variability 



Caswell et all (|l995bl ) showed that 12.2 GHz 



methanol masers commonly vary on time scales 
of three months by factors greater than 10 %, but 
rarely at the extreme levels that are seen in some 
water masers. On comparing the i r observations 



with th ose of ICaswell et al.1 (|19931 ). I Caswell et al 



(jl995bf ) noted that variations seen over several 3 
month periods were simi lar to the vari a tions o ver 
4 years. This prompted ICaswell et al.l (|1995a ) to 



suggest that the sources may go through quasi- 
periodic fluctuations. Indeed, long term moni- 
toring of some 12.2 GHz methanol masers have 
shown re markable periodicity in i ntens ity fluctua- 



tions (e.g. lGoedhart et al 



2005bll2009h. Observa 



tions of G 9.62+0.02 lead iGoedhart et al.l (|2005af ) 
to conclude, that at least in the case of this source, 
that the periodic flares were likely caused by peri- 
odic changes to the pump photon levels. 

We carried out our 12.2 GHz observations over 
two epochs separated by 6 months, observing a 
large fraction of sources at both epochs. However, 
repeat observations of sources were preferentially 
targeted towards those 6.7 GHz sources where no 
12.2 GHz methanol maser was observed in the first 
epoch, or those where only weak 12.2 GHz emis- 
sion was detected. The incomplete nature of the 
sample of sources studied at both epochs prevents 
an extensive study of variability from being com- 
pleted. Instead, sources with two epochs of data 
have been used to investigate the completeness of 
our search, as well as allowing a preliminary as- 
sessment of 12.2 GHz source variability. 

While approximately half of the 6.7 GHz 
methanol maser targets were observed during both 
epochs (264 of 580), only 56 of these resulted in 
a detection during either epochs. 28 of these 56 
sources could only be identified at one epoch, the 
strongest of which was 1.4 Jy when detected. The 



% of 6.7 GHz velocity range where the 6.7 GHz peak lies 



% of 6.7 GHz velocity range where the f 2.2 GHz peak lies 

Fig. 2. — Percentage of the 6.7 GHz methanol 
maser velocity range where the 6.7 GHz peak lies 
(top), and the 12.2 GHz lies (bottom). Only 
sources with 6.7 GHz methanol maser velocity 
ranges of more than 2 km s _1 are plotted. Here 
% corresponds to the minimum velocity and 100 
% corresponds to the maximum velocity. 
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average flux density of these once-detected-sources 
was ~0.6 Jy, which is comparable to the 5-u de- 
tection limit of our observations. However, as 
remarked in Section [31 many detections were con- 
fidently identified at the 3-cr level. In addition to 
this, averaging the observations of the 28 sources 
over the two epochs resulted in a decreased signal- 
to-noise in the resultant spectra. This shows that 
these 28 weak sources detected at only one epoch, 
do exhibit real variations in peak flux density of 
close to 50 % over the 6 month period: suggest- 
ing that the majority of these weak sources are 
more variable on average than the total 12.2 GHz 
methanol maser population. The luminosities of 
the 6.7 GHz methanol masers (discussed in Sec- 
tion I4.6.2[) that are the counterparts of these 28 
once-detected 12.2 GHz sources have an average 
luminosity that is almost five times lower than 
that of the full sample of 6.7 GHz masers with 
12.2 GHz counterparts. 

Given that half of the 56 sources that were ob- 
served during both epochs were detectable at only 
one epoch, it is likely that a significant number 
of additional 12.2 GHz methanol masers would be 
detected if a second epoch of observations were 
carried out for all sources. Of the 580 6.7 GHz 
methanol masers searched for 12.2 GHz methanol 
maser emission, 330 failed to yield detections and 
122 of these were observed only once. Since 12 
% (28 of 236) of the sources that were observed 
at both epochs but only revealed emission on one 
occasion, we can estimate that ~15 additional 
sources (12 % of 122) would be detected if a second 
epoch of observations carried out towards the 122 
sources that were only observed once and revealed 
no emission. 

We fail to detect three 12.2 GHz methanol 
masers that were previously detected by I Caswell et al 



(|!995bh (G 284.352-0.419, G 327.291-0.578 and 
G 349.067-0.017). These sources were reported to 
have peak flux densities of less than 1 Jy at the 
time of their detection. This is consistent with our 
estimation of the number of sources that would fall 
below our detection limits at any one epoch. 

Of the 28 sources that were observed and de- 
tected at both epochs, the average variation in 
peak flux density was 20 %. Only 8 of these 
sources had flux densities greater than 1 Jy at 
either epoch. Two sources showed variations of 
considerably more than 20 %, at levels of 71 and 



59 % respectively. Ten of the 28 sources showed 
no variation over the two epochs. 

4.5. Association with OH masers 

The locations of the target 6.7 GHz methanol 
masers (with and without 12.2 GHz counterparts) 
have been compa red with the lo cations of OH 
masers detected bv lCaswelll (|l998l ). ICaswelll (|l998l) 
conducted a complete search of a portion of the 
Galactic plane, along with a range of targeted ob- 
servations for OH masers. This has allowed us 
to not only identify methanol masers with associ- 
ated OH masers but also, for a large number of 
sources, those without. 393 of our 6.7 GHz maser 
target sample fall within the regions searched for 
OH masers by I Caswell (1998) and 144 of these are 
located within 1.5 arcsec of an OH maser. Using 
1.5 arcsec as an association threshold takes into 
account the positional uncertainties for both sets 
of observations and therefore allows the majority 
of methanol-OH associations to be identified while 
minimizing the number of chance associations. 
Furthermore, given that OH maser spots excited 
by a single object are often spread over region s 
greater than one arcsec ( Forster fc Caswell 1989), 
setting a tighter association threshold would not 
be appropriate. 



Both Lewis (2008) and Breen et all (|2010ah 



found that 12.2 GHz methanol masers were pref- 
erentially detected towards 6.7 GHz methanol 
masers that are associated with OH masers. 
Given that these searches were potentially bi- 
ased towards t he more evolved 6. 7 GHz methanol 
maser sources (|Breen et al.ll2010at ). we have tested 
this theory by comparing our complete methanol 
m aser sample wit h the complete OH maser sample 
of lCaswelll (|l998f) . 

We find, from our 12.2 GHz search towards our 
unbiased sample of 6.7 GHz methanol masers, that 
of the 393 6.7 GHz methanol masers with avail- 
able OH data, 144 have associated OH masers 
and that 86 of these also have detections at the 
12.2 GHz transition of methanol. Thus the search 
for 12.2 GHz towards 6.7 GHz methanol masers 
with OH masers has resulted in a detection rate 
of 60 %. For the 249 6.7 GHz methanol maser 
sources with no associated OH maser, 99 have as- 
sociated 12.2 GHz emission and therefore a search 
toward these 6.7 GHz sources has yielded a de- 
tection rate of 40 %. An alternative statistic is 
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that 46.5 % of the 185 6.7 GHz methanol masers 
with 12.2 GHz counterparts (that fall within the 
regions completely searched for OH maser emis- 
sion) have an OH maser counterpart. In compar- 
ison, the rate of OH maser association with the 
28 weak 12.2 GHz methanol masers detected at 
only one epoch (mentioned in Section f4.4p is much 
lower with only one association out of 12 sources 
that fall withi n the region completely searched by 
Caswell (1998). This finding is consistent with 
these weak variable sources representing perhaps 
the youngest 12.2 GHz sources. 

In the regions covered by complete searches for 
both 6.7 GHz methanol masers and OH masers 
( Caswell 19981) . the number of methanol masers 
detections out-number OH detections by more 
than two to one. Many estimates suggest that the 
association rate for meth anol masers to wards OH 
masers may surpass 80 % (|Caswel]ll2009h . However 
it is clear from comparison of the regions covered 
by both the OH maser survey of ICaswelll dl998h 
and the MMB catalog that the association rate 
for OH masers towards methanol masers is much 
lower, and in fact, is probably less than 40 %. 
These number statistics may imply that the life- 
time of methanol masers is much greater than that 
of the OH maser sources, but also shows that while 
the detection rate of 12.2 GHz methanol masers is 
greater at sites of OH masers, a search exclusively 
towards these sites would only uncover a fraction 
of the total population, since so many methanol 
masers are devoid of OH counterparts. The el- 
evated 12.2 GHz detection rate towards sources 
associated with OH masers also lends support to 
the idea that the 12.2 GHz sources are present to- 
wards the later stages of the period of evolution 
where 6.7 GHz masers are seen. 

4.6. Source luminosities 

Comparison between the luminosities of the 
two methanol maser transitions, and the other 
characteristics of each of the regions, has the 
potential to reveal much more insightful results 
than previous comp arisons of only flux densities. 
Breen et al.l (|2010a[) identified a relationship be- 
tween increasing 6.7 GHz methanol maser peak 
luminosity and evolution, whereby the more lu- 
minous 6.7 GHz methanol masers are associated 
with a generally later evolutionary stage of high- 
mass star formation. Comparison between am- 



monia linewidths and methanol maser luminosi- 
ties bvlWu et al.l (|2010l ) show that the lower lumi- 
nosity methanol masers correspond to the younger 
sources which is in agreement with the findings of 
i Breen et al l d . 2010ah . The analysis carried out by 
IBreen et al.l (|2010al ) was based on the peak flux 
densities of the 6.7 GHz sources extracted from 
the literature rather than the more appropriate 
integrated flux density, which is expected to be a 
more robust estimate of a larger volume of space. 
Furthermore, in many cases the observations of 
the two methanol maser transitions considered in 



Breen et al.l (|2010al ) were separated by 10+ years. 



Since we have integrated flux densities for both of 
the methanol maser transitions available for our 
large sample of sources, here we repeat som e of 
the analysis presented in lBreen et al.l (|2010al ) and 
offer new insights from additional investigations. 

4-6.1. Distances 

The calculated luminosity of the methanol 
masers are highly dependent on the adopted dis- 
tance measurements. We have calculated the kine- 
matic distances using the 6.7 GHz methanol maser 
peak velocities for our sample, using both t he re- 



vised prescription offered by iReid et al. (l2009h 



(which uses a flat rotation curve and updated 
solar parameters), as well as tha t of m ore tra- 
ditional model of brand fc Blit3 (|l993l ). The 
agreement between the peak methanol maser ve- 
locity and molecular line data is within a few 
km s 1 for most 6.7 GHz sources and are therefore 
adequat e for the computation of kinematic dis- 
tances (|Szvmczak. Bartkiewicz fc Richard's] 2007 : 



Pandian et al.ll2009f ) 

Comparison of the distances calculated using 
the two different methods shows that they are in 
gen eral very similar, a lthough those calculated us- 
ing Reid et al. (2003) tend to be slightly smaller. 
As such, the choice of model does not affect the 
calculated luminosities greatly. There are a few 
instances where the derived kinematic distances 
are not plausible, and in thes e cases neither the 
model of Brand & Blitd (|l993h or the prescription 
of Reid et al. (120091) are able to give an accept- 
able distance, indicating that either these sources 
have anomalous velocities or they lie in poorly con- 
strained regions of the Galaxy (especially when R 
< 3.5 kpc, where extreme departure s from circular 
rotation are known to be present (|Caswell et al 
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20101) 1. 

We adopt the prescription of Reid et al. ( 20091 ) 
which achieves better correspondence with par- 
allax measurements. Unfortunately sources with 
kinematic distance ambiguities have a more signif- 
icant impact on the luminosity distribution. The 
accurate assignment to the near or far distance re- 
quires techniques such as parallax measurements 
or Hi self absorption, which are beyond the scope 
of this paper. For ^15 % of sources the calculated 
distance is unambiguous, for those with near/far 
ambiguities we have used the near distances in all 
analysis. Methanol maser sources with calculated 
distances of less than 0.5 kpc or more than 20 kpc 
have been removed from this analysis as they are 
likely to be erroneous. Where we believe that out 
assumption of the near distance have affected our 
results we mention it in the text. 

4-6.2. Comparison between 6.7 GHz luminosities 
and Hi number densities 



Breen et al.1 (|2010al ) found that there was a re- 
lationship between the H 2 number densities of 
1.2 mm dust clumps (Hill et al. 2005) and the 
presence of UCH11 regions, with the more evolved 
1.2 mm dust clumps (i.e. those with associated 
UCH11 regions) showing lower densities than the 
younger dust clumps. Plots of 6.7 GHz peak 
luminosity versus H 2 number density (fig. 2 of 



Breen et al. (1201031) ) showed that the more lu- 
minous 6.7 GHz methanol masers were associ- 
ated with the less dense 1.2 mm dust clumps, im- 
plying a relationship between 6.7 GHz methanol 
source lu minosity and evoluti on. Furthermore, as 



noted in Breen et al 



(|2010ah . the appearance of 
the more evolved sources as being less dense may 
have been a conse quence of t he co nstant temper- 



ature assumed bv iHill et al . (2005), and may be, 
at least in part, indicating that the more evolved 
sources are warmer. In their analysis. iBreen et al 



(|2010al ) calculated the 6.7 GHz methanol maser lu- 
minosity using the peak flux density of the sources 
since the integrated flux densities were not avail- 
able in the literature. Here we use the luminosi- 
ties calculated from the integrated flux density to 
test whether the simpler analysis using luminosi- 
ties calculated from the peak flux densities led to 
any bias. F ig. [3] shows a similar plot to fig. 2 of 
but uses the integrated flux 



-20° (fr om the MMB ob servations) that overlap 
with the lHill et al. ( 2005 ) 1.2 mm dust clump sam- 
ple. Since the H2 number densities were calculated 
using kinematic distances for each source as pre- 
sented in Hill et al. (2005), (we have used the near 
distances in the case of ambiguities) to calculate 
the integrated 6.7 GHz luminosities for this plot 
for consistency. 

Linear regression analysis was carried out on 
the data contained in Fig. [3] and showed that the 
negative slope of the line of best fit of the distribu- 
tion is statistically significant, with a weak corre- 
lation between individual data points (correlation 
coefficient of 0.39). The null hypothesis that the 
line of best fit has a slope of zero can safely be 
rejected (p-value 0.007). The line of best fit has 
the equation 



log(L 6 . 7 ) = -0.72[0.25](log(H 2 )) + 5.52[1.17], 

where 'L6.7' is the integrated luminosity of 
the 6.7 GHz methanol maser (Jy km s -1 pc 2 ) 
and 'H2' is the H 2 number density (cm -3 ) cal- 
culated from the associated 1.2 mm dust clump 
emission. The numbers contained in the square 
brackets are the standard errors of the slope and 
the y-intercept res pectively. As mentioned in 
Breen et al. ( 2010al ). it is likely that at least some 
of the scatter present in the distribution of data 
points can be explained by the use of the near 
kinematic distances, which is an assumption prone 
to introducing errors. Additionally, since the 



1.2 mm dust continuum observations of Hill et al 



Breen et al 



densities of 6.7 GHz sources south of declination 



(2005) used to determine the H 2 number densities 
are of relatively low resolution, we are compar- 
ing the maser positions with the broader proper- 
ties of the associated region, as averaged by the 
large beam. A further contributor is likely to be 
the variable nature of the maser sources (see e.g. 
Goedhart. Gavlard fc van der Wald \2004 i. The 
correlation between data points is slightly lower 
for sources shown in Fig. [3] than those in fig. 2 of 
Breen et al. ( 2010al ). However, given the smaller 
sample size in Fig. [3] (46 compared to 113) it is dif- 
ficult to draw meaningful conclusions about what 
this may mean. The much smaller sample size 
is a result of th e fact that a large portion of the 
Hill et all |2005l) sampl e (towar ds which the obser- 
vations of IBreen et al.l (|2010al ) were targeted) lie 
north of the declination cut-off imposed on the tar- 
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get sources. Regardle ss, this analy s is show s that 



the result derived in iBreen et al.l (|2010al ) holds 



when the 6.7 GHz integrated luminosities are used, 
confirming that, in general, 6.7 GHz methanol 
masers increase in luminosity as they evolve. 



4-6.3. The relationship between the 
12.2 GHz source luminosities 



6. 7 and 



Breen et al.l (|2Q10al ) identified a relationship 
between the ratio of 6.7 to 12.2 GHz peak flux den- 
sities with the luminosities of the respective data. 
The ratio of 6.7 to 12.2 GHz peak flux density ap- 
parently increased with increasing 6.7 GHz peak 
luminosity, implying that as the source evolved, 
the 6.7 GHz methanol maser increased in flux den- 
sity more rapidly than the associated 12.2 GHz 
methanol maser. We have repeated that analysis 
on our much larger and complete sample using not 
only the peak flux densities of the sources, but also 
comparing the integrated flux densities. 

Log-log plots of: 

1. the ratio of 6.7 GHz to 12.2 GHz peak flux 
density ratio versus the 6.7 GHz integrated 
luminosity; 

2. the ratio of 6.7 GHz to 12.2 GHz methanol 
maser integrated flux densities versus the 
6.7 GHz integrated luminosity; 

3. the ratio of 6.7 GHz to 12.2 GHz methanol 
maser peak flux densities versus the 12.2 GHz 
integrated luminosity, 

were created and inspected. Together with linear 
regression analysis, these plots showed simila r re- 
sults to those reported in IBreen et al.l (|2010al ). In 
all cases the slopes of the line-of-best-fit were con- 
sistent with the notion that the luminosity of both 
of the methanol maser transitions increase with 
evolution and that the rate of increase is lower for 
the 12.2 GHz sources. 

Table[2]shows a summary of the results from the 
linear regression analysis for the data correspond- 
ing to the three plots described above. From this 
table it can be seen that the correlation coefficients 
in all cases correspond to a weak correlation be- 
tween individual data points and this scatter can 
largely be accounted for by the inaccuracies intro- 
duced by the assumption that all of the methanol 
masers lie at the near kinematic distance. 



"a 

S 
to 



CD 

CD 









o 




o 








o 








o 






o 










° I 












oo 


o 




( 


3 o o o 




o °o 










oo 






o 




oo ^\ 






o 



















0° 























o 


o 



3.5 4.0 4.5 5.0 5.5 



Log H_2 number density 



Fig. 3. — Log of the 6.7 GHz integrated luminosity 
(Jy km s _1 kpc 2 ) versus the log of the H2 num- 
ber density of the associated 1.2 mm dust clump 
(cm -3 ). The line of best fit is also plotted. 



Table 2: Relationship between relative 6.7 and 
12.2 GHz luminosities, searching for any depen- 
dence on 6.7 or 12.2 GHz integrated luminosities, 
and described in the text as plots 1, 2 and 3. 



Plot 


Slope 


Error 


p- value 


R 


1. 


0.11 


0.03 


2.28 x e" 4 


0.23 


2. 


0.11 


0.03 


1.79 x c~ 3 


0.19 


3. 


-0.12 


0.03 


6.24 x e" 5 


0.25 



Column 2 shows the slope of the line of best fit, 
column 3 shows the standard error associated with 
the slope of the line, column 4 shows the p- value 
relating to the null hypothesis that the line has 
a zero slope, and column 5 gives the correlation 
coefficient. 
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Log 6.7 GHz integrated luminosity 



^ i i i r 

-2 2 4 6 

Log 6.7 GHz integrated luminosity 

Fig. 4. — Log 12.2 GHz integrated luminosity ver- 
sus the log of the 6.7 GHz integrated luminosities 
(luminosities in units of Jy km s _1 kpc 2 ). Over- 
laid is the line of best fit (dashed) and line with a 
slope of 1 (solid). 



Fig. 5. — Log 12.2 GHz integrated luminosity ver- 
sus the log of the 6.7 GHz integrated luminosities 
(luminosities in units of Jy km s _1 kpc 2 ). The 
unfilled circles represent those sources with an as- 
sociated OH maser and the filled circles show the 
sources with no OH maser counterpart. Overlaid 
is the line of best fit for sources with OH masers 
(dashed), sources without OH masers (dotted), 
the full sample with OH information (dot-dashed) 
and line with a slope of 1 (solid) for comparison. A 
summary of the linear regression analysis for each 
of the source categories is shown in Table [3] 
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A comparison between the 6.7 GHz and 12.2 GHz 
luminosities calculated from their integrated flux 
densities is shown in Fig. |4] Not surprisingly, this 
figure shows that the majority of these methanol 
maser sources have intermediate luminosities, ev- 
ident by the concentration of sources in the log 
luminosity range 1 to 3 and to 2 for the 6.7 GHz 
and 12.2 GHz sources respectively. The relation- 
ship between the 6.7 and 12.2 GHz methanol 
maser luminosities is well estimated by the line 
of best fit as it shows a relatively tight corre- 
lation (correlation coefficient is 0.85) between 
individual data points, with no extreme out- 
liers. This plot similarly shows that for higher 
6.7 GHz luminosities, the ratio of 6.7 GHz to 
12.2 GHz source flux density increases. This is 
consistent with the maser mode ls presented in 
Cragg. Sobolev fc Godfrey! (j2005h wherein fig. 2 
shows that the 6.7:12.2 GHz flux density increases 
as the source evolves. The line of best fit has the 
equation 

log(L 12 . 2 ) = 0.89[0.03](log(L 6 . 7 )) - 0.90[0.09], 

where 'L12.2' is the integrated luminosity of the 
12.2 GHz methanol maser and 'L6.7' is the inte- 
grated luminosity of the 6.7 GHz methanol maser. 
The numbers contained in the square brackets are 
the standard errors of the slope and the y intercept 
respectively. 

We have further investigated this relationship 
by splitting the data into two categories. The first 
category includes all 6.7 and 12.2 GHz methano l 
masers with associated OH masers (|Caswelllll998l ) 
and the second includes those sites with both 6.7 
and 12.2 GHz methanol masers but no associ- 
ated OH maser. Fig. [5] presents a plot similar 
to Fig. S] but with methanol sources with asso- 
ciated OH masers distinguished from those with- 
out. Methanol sources that have not been sensi- 
tively searched for OH masers have been excluded 
from the plot (a total of 58 sources). Interestingly, 
the methanol sources with OH maser counterparts 
tend towards the higher luminosities, lending fur- 
ther support to the idea that methanol masers 
tend to increase in luminosity as they evolve (since 
OH masers are known to be associated with a 
later stage of evolution ( Forster fc Caswell 1989t 
Casweltifl99"7K 

The summary of the linear regression analysis 



carried out on the three categories of sources dis- 
played in Fig. [5] (all sources, sources with OH and 
sources with no OH counterpart) are presented 
in Table [3] The characteristics computed in the 
linear regression analysis for all sources shown in 
Fig. [5] are almost identical to those derived for the 
complete sample (Fig. [4]). Sources with associated 
OH masers seem to have a slightly more positive 
slope (~1) than the sources with no OH maser 
counterparts that have a slope equal to that of the 
full sample (0.90). Taking into account the stan- 
dard errors associated with these two lines shows 
that they could have the same slope (i.e. both 
less than one), although the possibility that this 
is a real difference is intriguing. If the difference is 
real, it implies that although the methanol maser 
luminosity increases more rapidly in the 6.7 GHz 
tranisition as the source evolves, towards the evo- 
lutionary stage where OH masers are present this 
effect lessens. 

4-6. 4- Preliminary analysis of the luminosity of 
individual maser features 

A preliminary analysis has been carried out on 
individual features of 6.7 GHz masers with respect 
to their 12.2 GHz counterparts. Since the two 
methanol maser transitions are t ypically co-spatial 
to w ithin a few miliarcseconds (jMoscadelli et al. 
20021) . and that this is especially likely when the 



spectral profiles of the res pective transitions are 
similar (jNorris et al.lll993l ). we are able to mean- 
ingfully draw these comparisons even in the ab- 
sence of detailed high-resolution maps. Future 



Table 3: Summary of results from linear regression 
analysis for data presented in Fig. [5j 



Category 


Slope 


Error 


p- value 


R 


full sample 


0.90 


0.05 


<2 x e" ie 


0.82 


with OH 


1.01 


0.08 


<2 x e" 16 


0.80 


no OH 


0.90 


0.06 


<2 x e" 16 


0.83 



Column 1 gives the category of the sources in the 
sample, column 2 shows the slope of the line of 
best fit, column 3 shows the standard error asso- 
ciated with the slope of the line, column 4 shows 
the p- value relating to the null hypothesis that the 
line has a zero slope, and column 5 gives the cor- 
relation coefficient. 
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studies of these sources will include an analysis 
of the individual features of a much large group of 
sources. 

Analysis of the individual features associated 
with each maser source has the potential to re- 
veal much more about the sources, as well as 
the individual features, which show emission at 
both methanol transitions, or are devoid of a 
12.2 GHz counterpart. Since the age of the excit- 
ing star will be common to individual maser fea- 
tures, their analysis will reap greater information 
on the study of the relative evolutionary stage that 
each methanol transition is associated with. Fur- 
thermore, insights into the conditions favourable 
to both transitions can be achieved by carrying 
out the analysis on such a large sample of sources. 

Fig. |6] displays the peak luminosity of each 
feature associated with five 6.7 GHz methanol 
masers, versus the ratio of 6.7 to 12.2 GHz lu- 
minosity of the corresponding feature. The five 
sources presented in this plot all exhibit three or 
more 12.2 GHz spectral features and were cho- 
sen randomly from this population. In addition 
to the ratios calculated from detections at both 
transitions (in individual features), upper limits 
on the ratios of 6.7 GHz features with no 12.2 GHz 
counterpart, are shown. The individual features of 
the five sources are colour-coded, revealing several 
trends. 

Firstly, the 6.7 GHz luminosity versus the ratio 
of 6.7 to 12.2 GHz luminosity of each spectral fea- 
ture are clustered together in groups containing all 
the features associated with the same source. It is 
expected from Sect ion 14.6.21 fin conjun ction with 
results presented in lBreen et al. I (|2010ah l that the 
ratio of 6.7 to 12.2 GHz luminosity will increase 
with 6.7 GHz luminosity, and therefore evolution 
of the source. The clustering of the features as- 
sociated with each source supports this idea since 
the exciting source will be the same, and therefore 
at the same evolutionary stage, for all spectral fea- 
tures. 

A second implication from Fig. |6] comes from 
the location of the lower limits for the ratios. It 
can be seen that, while the features which show 
emission at both 6.7 and 12.2 GHz are clustered, 
the upper limits for features with no detectable 
12.2 GHz counterpart lie on the outskirts of the 
cluster; since these are upper limits, if anything, 
they are likely to lie even further from the cluster. 
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Fig. 6.— Log of the 6.7:12.2 GHz individual fea- 
ture ratio versus the log 6.7 GHz peak luminos- 
ity for every spectral feature associated with five 
methanol maser sources. Dots represent 6.7 GHz 
features with a 12.2 GHz counterpart, while tri- 
angles show features with no detectable 12.2 GHz 
counterpart, thus showing the lower limit of the 
6.7:12.2 GHz feature ratio. Colours distinguish 
between individual maser sites. The dashed line 
shows the line of b est fit from the sa mple of 
sources presented in iBreen et al. I (l2010ah . fig. 2, 
where a plot of 6.7 GHz peak luminosity versus 
6.7:12.2 GHz peak ratio is shown for the peak 
6.7 GHz feature associated with each source. 
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This implies that many of the features with no 
12.2 GHz counterpart that are present on this plot 
are unlikely to be just below the detection limits 
of the search. 

Also shown in Fig.[6]is the line of best fit calcu- 
lated from the linear regression carried out on the 
lar ger group o f source s shown on a similar plot 
in iBreen et al. ( 2010af ). fig. 2 which shows only 
the peak feature for each of the methanol maser 
sources. The plotted line is the same (to within 
the errors) as the line of best fit calculated using 
the points plotted on Fig. |6] (omitting the upper 
limits) . 

4- 6. 5. Methanol maser luminosity and velocity 
range comparison 

There is a striking difference in the values 
of both luminosity and velocity range of the 
methanol masers in the association categories: 
firstly, the 6.7 GHz sources devoid of both 12.2 
and OH maser counterparts; secondly, those with 
associated 12.2 GHz emission but not OH; and 
thirdly those with both 12.2 GHz and OH maser 
counterparts. We find that the 12.2 GHz methanol 
maser sources with associated OH masers have 
higher luminosities and velocity ranges than those 
with no OH counterpart. OH maser associations 
(or lack thereof) with the methanol masers have 
been determined from comparisons w ith the OH 
maser sources listed in Cas well (Il998h . 

Table [4] shows the average and median veloc- 
ity ranges and luminosities of both the 6.7 and 
12.2 GHz methanol masers broken up in cate- 
gories according to their associations. In the 
top half of the table the characteristics of the 
6.7 GHz methanol masers are given and show 
a progressive widening of velocity range and in- 
creasing luminosity from the category which in- 
cludes 6.7 GHz methanol masers with no 12.2 GHz 
methanol or OH maser counterpart; to 6.7 GHz 
methanol masers with associated OH maser emis- 
sion but no 12.2 GHz methanol maser emission; to 
6.7 GHz methanol masers with associated emission 
at the 12.2 GHz transition but no associated OH 
maser; and finally those 6.7 GHz methanol masers 
with both 12.2 and OH maser counterparts. The 
bottom half of the table gives the characteristics 
of the 12.2 GHz methanol masers with and with- 
out associated OH masers. In this case, a sim- 
ilar trend is seen, with the sources without OH 



masers showing smaller 12.2 GHz maser velocity 
ranges and lower luminosities than those 12.2 GHz 
sources with associated OH masers. 

Figs [7] and |5] show plots of log 6.7 GHz luminos- 
ity versus 6.7 GHz velocity range and log 12.2 GHz 
luminosity versus 12.2 GHz velocity range, respec- 
tively. In these figures, there is one plot of the 
maser luminosity calculated using the integrated 
flux density and one using the maser peak flux 
densities. Fig. [7] shows that there is an overall 
trend between increasing 6.7 GHz methanol maser 
luminosity (for both integrated and peak) and ve- 
locity range, with a significantly positive slope 
and a moderate correlation between data points. 
These plot also highlight the segregation of the 
values correspon ding to the d i fferent association 
categories. Like Breen et al. ( 2010al ). we argue 
that, in general, more evolved 6.7 GHz methanol 
masers exhibit higher luminosities. The values re- 
ported in Table 2] along with the data shown in 
Fig. [7] support this idea. This plot offers addi- 
tional insight, with a clear tendency for the more 
evolved 6.7 GHz methanol masers to show larger 
velocity ranges. The figure showing the lumi- 
nosities calculated from the integrated flux den- 
sities suggests that the gas volume conducive to 
6.7 GHz methanol maser emission also increases 
as the source evolves. Alternatively, the increase 
in velocity range may, at least partially, be at- 
tributed to increased internal motions within the 
more evolved sources. 

There is, of course, the possibility that this 
relationship is influenced by sensitivity, i.e. for 
example, that the 6.7 GHz sources with associ- 
ated 12.2 GHz methanol and OH masers have 
higher flux densities, increasing the likelihood that 
a more extensive velocity range would be iden- 
tified. However, Fig. [6] suggests that sensitivity 
biases are not the culprit for the difference be- 
tween the sources. Even if this is the case for 
some sources, Fig. [5] shows that the characteris- 
tics of the 6.7 GHz methanol masers change with 
evolution regardless (since this shows that the flux 
density increases). Furthermore, if sensitivity lim- 
itations were precluding us from observing the full 
extent of the 6.7 GHz methanol maser velocity 
range, then we would expect this to affect all of 
the association categories equally since the differ- 
ent categories will be interspersed through a range 
of distances in the Galaxy. A sample of the outly- 
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Table 4: Average and median velocity ranges and luminosities of 6.7 and 12.2 GHz methanol masers broken 
into categories on the basis of their methanol and OH maser associations. 



Associations 


Sample 


Velocity range 


Integrated luminosity 




size 


(km 




(Jy km s 


1 kpc 2 ) 






average 


median 


average 


median 


6.7 GHz 












no 12.2 or OH 


136 


5.2 


3.9 


70 


25 


no 12.2 but OH 


56 


6.1 


5.5 


269 


64 


12.2 but no OH 


93 


8.3 


7.6 


465 


142 


12.2 and OH 


84 


10.0 


9.4 


3819 


892 


12.2 GHz 












no OH 


93 


3.1 


1.6 


35 


17 


OH 


84 


4.7 


3.0 


667 


73 



In the top half of the table, characteristics of the 6.7 GHz methanol masers are shown, broken into their 
association categories, while the second half of the table presents the characteristics of the 12.2 GHz methanol 
masers. Column 1 specifies the association category for the 6.7 (top) and 12.2 GHz (bottom) methanol 
masers; column 2 gives the number of sources in each category; column 3 gives the average velocity range 
(km s _1 ); column 4 gives the median velocity range (km s _1 ); column 5 gives the average luminosity (Jy 
km s _1 kpc 2 ); and column 6 gives the median luminosity (Jy km s _1 kpc 2 ). 



ing sources shown in top plot of Fig. [7] have been 
investigated, specifically those sources represented 
by purple dots and black triangles centred on ^1 
on the x-axis and ~17 on the y-axis. These sources 
show unexpectedly large velocity ranges for their 
calculated luminosities. An investigation of these 
sources shows that they all have low flux densities 
(<4 Jy) and this together with their location on 
the plot implies that they might actually lie at the 
far distance rather than the near distance that has 
been assumed. Therefore, a comparison between 
the integrated flux densities and velocity ranges 
of sources may provide a viable, indirect way of 
resolving near/far distance ambiguities. 

A further group of sources are presented in 
Fig. [3 and Table HJ these are those 6.7 GHz 
methanol masers with associated OH maser emis- 
sion but no 12.2 GHz methanol maser emission 
(represented by green asterisks). These sources 
do not fit in with the evolutionary sequence that 
seems to describe the majority of our sources. The 
6.7 GHz methanol sources in this category have av- 
erage and median velocity ranges and luminosities 
that fall in between the categories of 'no 12.2 or 
OH' and '12.2 but no OH'. It seems, therefore, that 
for some sources, the OH masers switch on prior to 
the 12.2 GHz methanol masers. It can be seen that 
these sources basically populate the same range as 



that of the two categories of 'no 12.2 or OH' and 
'12.2 but no OH'. These sources need further inves- 
tigation to understand their nature. Of particular 
interest would be to compare the abundances of 
methanol in these sources compared to the other 
categories, as well as determining their associa- 
tions with UCHu and hyper-compact Hn regions. 

Fig.[8]displays the characteristics of the 12.2 GHz 
methanol masers, broken into the categories of 
sources associated with an OH maser (red crosses) 
and sources with no OH maser counterpart (black 
triangles). These plots similarly show evidence 
that the 12.2 GHz methanol masers increase in 
velocity range and luminosity as they evolve, since 
OH masers are f ound at a later stage of h i gh-mass 
star formation (jForster fc Caswelll Il989t ICaswelll 
Il997j ). The congregation of a number of sources 
with small velocity ranges and the log of the inte- 
grated luminosity between and 2, perhaps shows 
that the 12.2 GHz sources are more affected by 
sensitivity limitations prohibiting the full velocity 
ranges to be identified. The sensitivity limits of 
the 6.7 and 12.2 GHz searches were similar, and 
together with the expectation that the associated 
12.2 GHz methanol masers sources are weaker, 
means that this is not unexpected. However, anal- 
ysis of individual features associated with five 6.7 
and 12.2 GHz methanol masers shown in Sec- 
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Fig. 7. — Log of the 6.7 GHz integrated luminosity (Jy km s -1 kpc 2 ) versus the 6.7 GHz methanol maser 
velocity range (km s _1 ) (top) and log of the 6.7 GHz peak luminosity (Jy kpc 2 ) versus the 6.7 GHz methanol 
maser velocity range (km s _1 ) (bottom). To be present on these plot sources must have acceptably deter- 
mined kinematic distances and integrated/peak flux densities. Shown on the plots are four categories of 
sources: 6.7 GHz sources with no associated 12.2 GHz methanol or OH masers (136 purple dots); 6.7 GHz 
methanol masers with detectable 12.2 GHz methanol masers but no associated OH masers (93 black trian- 
gles); 6.7 GHz methanol masers with both 12.2 methanol and OH maser counterparts (81 red crosses); and 
the 6.7 GHz methanol masers with no 12.2 GHz methanol maser counterpart but with OH maser emission 
(56 green asterisks). Overlaid are the lines of best fit that corresponds to all of the 6.7 GHz masers presented 
on each plot (linear regression gives: (top plot) slope=2.57[0.27], p- value below 2 x e -16 , and a correlation 
coefficient of 0.48; (bottom) slope=2.03[0.29], p-value=6.31 x e -12 , and a correlation coefficient of 0.34). 
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Fig. 8. — Log of the 12.2 GHz integrated luminosity (Jy km s _1 kpc 2 ) versus the 12.2 GHz methanol maser 
velocity range (km s _1 ) (top) and log of the 12.2 GHz peak luminosity (Jy kpc 2 ) versus the 12.2 GHz 
methanol maser velocity range (km s ) (bottom). Shown on the plots are 12.2 GHz sources with no 
associated OH masers (93 black triangles); and 12.2 methanol masers with OH maser counterparts (84 red 
crosses). Overlaid are the lines of best fit that corresponds to all of the 12.2 GHz masers presented on each 
plot (linear regression gives: (top) slope=1.86[0.31], p-value=1.72 x e~ 8 , and a correlation coefficient of 0.40; 
(bottom) slope=1.82[0.33], p-value=1.40 x e -7 , and a correlation coefficient of 0.38). 
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tion 14.6.41 implies that, at least for these sources, 
the 6.7 GHz features with no detectable 12.2 GHz 
counterpart are unlikely to be a result of sensitiv- 
ity limitations. Both scenarios are probably true 
for some number of sources. 

From analysis of the ratio of integrated flux 
densities of the two methanol transitions in Sec- 
tion 14.31 it was suggested that there was a larger 
volume of gas that was conducive to 6.7 GHz 
methanol than their 12.2 GHz counterparts. Data 
presented in Table [4] and Figs [7] and [8] supports 
this, with a clear tendency for the 12.2 GHz 
methanol maser sources to have less extensive ve- 
locity ranges. In addition, this table and fig- 
ures show strong evidence for the gas volume to 
increase as the sources evolve (for both the 6.7 
and 12.2 GHz sources). In Fig. 01 and discus- 
sion thereof, we show that the luminosities of both 
the 6.7 GHz methanol maser and the 12.2 GHz 
methanol mascrs increased with evolution, but 
the the rate of increase was much greater for the 
6.7 GHz sources. The slopes of the lines derived 
from linear regression analysis of the data pre- 
sented in Figs [7] and [8] show that a similar effect 
is seen in the velocity ranges of the two methanol 
transitions. The slope associated with the 6.7 GHz 
methanol sources is 1.4 times that associated with 
the 12.2 GHz sources, showing that the velocity 
range of the 6.7 GHz methanol sources increases 
more rapidly than that of their 12.2 GHz methanol 
maser counterpart. 

These findings support the idea that the differ- 
ent combinations of 6.7 GHz methanol, 12.2 GHz 
methanol and OH masers, trace differing evo- 
lutionary stages of high mass star formation in 
the following sequence (starting from youngest); 
6.7 GHz methanol masers devoid of both 12.2 GHz 
counterparts and OH masers; 6.7 GHz methanol 
masers with 12.2 GHz counterparts but no asso- 
ciated OH maser emission; and finally, 6.7 GHz 
methanol masers with both 12.2 GHz counterparts 
and associated OH masers. This interpretation is 
self consistent with all other analyses completed 
and is also consistent with the independent re- 
~ (|2010ah . 



suits of Breen et al 



Further to this re- 
sult, we find that not only the luminosity of the 
6.7 GHz methanol maser increases more rapidly 
than that of the associated 12.2 GHz maser, but 
also the 6.7 GHz velocity range. This supports 
an evolutionary scenario whereby both the lumi- 



nosity and the gas volume conducive to the 6.7 
and 12.2 GHz masers increases as they evolve (al- 
though at a lesser rate for the 12.2 GHz sources). 

4.7. Summary of evolutionary results 

In this section we give a summary of the 
main evolutionary arguments presented in the 
previous sections and highlight the consistency 
of the these results. Our findings are consistent 
with those that are well established in the litera- 
ture; particularly that OH ma sers are associated 
with a later stage of evolution ( Forster fc Caswell 
19891: ICaswelllll997h. Also, as previously men- 



tioned. IWu et all (|2010t ) found that low luminos- 
ity 6.7 GHz methanol masers are associated with 



small ammonia linewidths, which Longmore et al 



(|2007l ) suggests are indicative of young source; 
consistent with the methanol masers increasing 
in flux den s ity wit h evolution as suggested by 



Breen et all ( 2010a ) and consistent with our cur- 



rent findings. 

Below we list some of integral findings and their 
evolutionary interpretations, followed by some of 
the main supporting evidence for each of the state- 
ments. 

12.2 GHz detection rate is higher towards 
the more evolved sources. This is evident from 
our comparison of detection rates towards sources 
with and without OH maser emission, but also 
in the high detection rates of previous searches 
which were biased towards more evolved star for- 
mation regions. This suggests that 12.2 GHz 
methanol masers are associated with second half 
of the 6.7 GHz methanol maser lifetime. This is 
further supported by the fact that there are few 
associations between OH mascrs and very weak 
12.2 GHz methanol masers (which we interpret as 
being the youngest 12.2 GHz sources; see the next 
point). 

The methanol maser emission increases 
in luminosity with evolution. Comparison 
between the Ho number densities calculated by 
iHill et all (|2005f ) a nd methanol maser luminosity 
(both here and in Breen et al ( 2010al )) shows a 
trend that is consistent with the sources increasing 
in luminosity with evolution. Comparison between 
the luminosity of the two methanol maser tran- 
sitions shows that the 6.7 GHz methanol maser 
emission is increasing at a higher rate than the 
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12.2 GHz sources. Furthermore, sources that show 
methanol maser have higher luminosities when 
they are also associated with an OH maser. 

The methanol maser emission have wider 
velocities ranges as they evolve. Luminosity 
of both methanol maser transitions increases with 
increasing velocity range, showing that the veloc- 
ity range of sources is also increasing with evolu- 
tion. Furthermore, both the velocity range and 
luminosity of sources increases in the association 
groups. In the case of the 6.7 GHz sources this 
increase is from solitary 6.7 GHz sources, to those 
with 12.2 GHz counterparts, and finally those with 
both 12.2 GHz counterparts and OH masers. 



Ther efore we suggest, similarly to iBreen et al 



(|2010al ) but with a few additions, that 6.7 GHz 
methanol masers devoid of 12.2 GHz or OH maser 
emission present the youngest sources, followed by 
those with associated 12.2 GHz masers and finally 
those sources which exhibit both methanol maser 
transitions as well as OH maser emission. We fur- 
ther suggest that the methanol masers increase 
in both luminosity and velocity range with evolu- 
tion and that this affect is greater in the 6.7 GHz 
sources. 

4.8. Association with GLIMPSE sources 

The locations of the detected 12.2 GHz methanol 
masers, along with the 6.7 GHz sources with no 
detectable 12.2 GHz counterpart have been com- 
pared to both the GLIMPSE point source cata- 
log and the catalog o f Extended Green Object s 



(EGOs) presented in ICveanowski et ah (2008) 



The characteristics of infrared sources have been 
regarded as an excellent determinant of relative 
source ages over relat ively long time-scales (e.g. 



Rathborne et alj |2010). The comparison between 



such a large number of methanol maser sources 
with the mid-infrared data is an important step 
in comparing the relative evolutionary stages as- 
sociated with different maser species. 

4-8.1. Association with GLIMPSE point sources 

520 of our target 6.7 GHz methanol fall within 
either the GLIMPSE I or GLIMPSE II regions. 
Using an association threshold of 2 arcsec, 316 
of our target methanol maser sources are coinci- 
dent with a GLIMPSE point source listed in the 
catalog (a detection rate of 61 %). The asso- 



ciation rate for 6.7 GHz methanol masers with- 
out associated 12.2 GHz emission is higher (192 
of 291, or 66 %), than for sites exhibiting the 
two methanol maser transitions (127 of 229, or 
55 %). Methanol maser sources coincident with 
OH masers have a similar association rate as for 
the sources showing emission at 12.2 GHz, a rate 
of 57 % (75 of 132). The similarity of the detection 
rates for sources with 12.2 GHz methanol masers 
or OH masers is expected given the high rate of 
coincidence for 12.2 GHz masers towards sites of 
6.7 GHz methanol and OH masers. 

Ellingsenl (|2006h found a GLIMPSE point 



source coincidence rate of 56 % (23 of 41) for a 
subset of their 6.7 GHz sources that were detected 
in a complete search, and this rate dropped to 
52 % (29 of 56) when the full sample was consid- 
ered. These association rate statistics, considered 
in conjunction with our own, suggests that the rel- 
atively more evolved star formation regions have 
fewer associations with GLIMPSE point sources. 
Evidence for this includes the fact that our sample 
of 6.7 GHz methanol masers yield the highest as- 
sociation rate (thereby including a high number of 
weak, possibly younger, sources) and this rate gets 
progressively less when stronger 6.7 GHz methanol 
masers or those associated with 12.2 GHz or OH 
maser emission, are considered. Since it is more 
likely for the more evolved sources to no longer 
be a point source in the GLIMPSE bands a lower 
association rate is expected. 

(l2010l) 



Gallawav et al 



have undertaken ex- 
tensive analysis of a near complete sample of 
methanol masers with associated GLIMPSE sources 
Ordinarily analyses such as these have been lim- 
ited to sources coincident with GLIMPSE point 
sources and to those IRAC bands with measured 
fluxes (which limits the sample). iGallawav et al 



( 2010l ) have used an adaptive non-circular aper- 
ture photometry technique to measure the fluxes 
of GLIMPSE sources towards 512 methanol 
masers. From this large sample of comparisons, 
they show that the r esults of the smaller sample 
studies bv lEllingsenl ( 2006 ) remain valid. 

Fig. [9] shows a sample of GLIMPSE colour- 
colour and colour-magnitude diagrams for the 
6.7 GHz methanol masers coincident wit h GLIMPS E 
point sources. We find, similarly to lEllingsen 
(120061 ). that the maser associated sources are sig- 
nificantly redder than the majority of the compar- 
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ison sources (all of the GLIMPSE point sources 
within 30 arcmin of I = 326.5°, b = 0°) on the 
plot, lying in significantly different regions in each 
case. Evident from the [5.8]-[8.0] versus [3.6]-[4.5] 
colour-colour plot, the methanol maser associated 
sources appear to have a large relative excess of 
4.5 /im flux. Interestingly, this plot along with 
the [3.6]-[4.5] versus 8.0 /im colour-magnitude di- 
agram show that the majority of the sources that 
have similar 3.6 and 4.5 /im magnitudes are with- 
out associated 12.2 GHz emission. Aside from 
this, there are few obvious differences between the 
GLIMPSE point sources associated with both 6.7 
and 12.2 GHz methanol maser emission and those 
associated with only 6.7 GHz methanol maser 
emission. 

Additional plots similar to those presented in 
Fig. [9] were constructed with the methanol masers 
sources associated with OH maser emission distin- 
guished from those with no associated OH maser 
emission. The sources with associated OH masers 
are relatively brighter in the 4.5 and 8.0 /im bands 
than those 6.7 G Hz metha n ol ma sers with no OH 
maser emission. lEllingsenl (|2006l ) similarly noted 
that 6.7 GHz methanol masers with OH maser 
counterparts tended to have brighter 8.0 /im mag- 
nitudes. 



Ellingsenl (|2006l ) commented that the majority 



of masers associated with GLIMPSE point sources 
meet the criteria [3. 6] -[4. 5] > 1.3 and 8.0 /im 
mag < 10 as well as [3.6]-[8.0] > 4.0 mag. From 
Fig. [9] it can be seen that these criteria are met by 
most of our methanol maser sources, although per- 
haps fail to capture some of the methanol maser 
sources, particularly a number of sources with no 
12.2 GHz maser counterpart. From our investiga- 
tion it seems that these criterion could be refined 
to [3.6]-[4.5] > 0.8, 8.0 /mi mag < 10 and [3.6]- 
[8.0] > 2.0 mag. This may imply that a number 
of the younger sources have a less excessive 4.5 to 
3.6 /tm flux and likewise for 8.0 /im band. 

Analysis of the GLIMPSE properties associ- 
at ed with the somewha t smaller sample presented 
in iBreen et al. ( 2010a ) resulted in the conclusion 
that the masers themselves are much more sensi- 
tive to evolutionary changes than the associated 
mid-infrared emission. Apart from the few trends 
mentioned above, we similarly draw this conclu- 
sion. 



4-8.2. Association with Extended Green Objects 
(EGOs) 

The locations of the 6.7 GHz methanol maser 
sources were compared with the positio ns of EGOs 
presented in Cvganowski et al. ( 20081) . We have 
adopted an association threshold of 10 arcsec to 
ensure that all associations were captured while 
minimising the number of chance associations. We 
find that 97 of the 580 6.7 GHz MMB methanol 
masers are associated with EGOs, and of these 
sources 49 have 12.2 GHz methanol counterpart. 
This equates to an overall 12.2 GHz detection rate 
of 50.5 % (49 of 97) for 6.7 GHz methanol masers 
associated with EGOs and 42.2 % (121 of 287) 
towards those without (after discarding 6.7 GHz 
sources that lie outside the regions inspected by 
Qyganowski et all (120081) for EGOs; 10° <1< 65° 
and 295° < I < 350°, b =±1°). This result means 
that in the case of a 6.7 GHz methanol maser 
that is associated with an EGO, the detection of a 
12.2 GHz methanol maser is more likely than to- 
wards a 6.7 GHz methanol maser that is not asso- 
ciated with an EGO. The 12.2 GHz detection rates 
towards EGOs is statistically significantly higher 
than the detection rate for the entire 6.7 GHz sam- 
ple which is 43 %. 

These comparisons also show that when tak- 
ing into account only those 6.7 GHz methanol 
masers that fall within the regions inspected by 
Cvganowski et all (|2008l) that 25.3 % (97 of 384) of 
6.7 GHz methanol masers have associated EGOs 
while this increases to 28.8 % (49 of 170) when 
you take into account only those 6.7 GHz masers 
that have associated 12 .2 GHz methanol masers. 
Cvganowski et al. I (l2009h found that the detection 
rate of 6.7 GHz methanol masers towards was 
much higher at a rate of ~64 % and suggested 
that this detection rate was higher than targeted 
searches that had used other selection criteria. 
While the detection rate is high, it is clear that 
a 6.7 GHz methanol maser search exclusively tar- 
geting EGOs may only recover approximately one- 
quarter of the total population. 

69 of the 6.7 GHz methanol masers that are as- 
sociated w ith EGO s have been searched for OH 
masers by I Caswell ( 1998 ). 39 of these sources 
have an associated OH maser, an association rate 
of 57 %. In the case where an EGO is asso- 
ciated with both a methanol and an OH maser 
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Fig. 9. — GLIMPSE colour-colour and colour-magnitude diagrams. Clockwise from top left: [5.8]-[8.0] 
versus [3.6]-[4.5] colour-colour diagram; [3.6]-[8.0] versus [3.6]-[5.8] colour-colour diagram; [3.6]-[4.5] versus 
8.0 /im colour-magnitude diagram; and [5.8]-[8.0] versus 5.8 /jm colour-magnitude diagram. GLIMPSE point 
sources associated with sources showing emission at both the 6.7 and 12.2 GHz methanol maser transitions 
are represented by blue squares and those associated with only the 6.7 GHz transition are represented by 
red circles. In each plot the black dots represent all of the GLIMPSE point sources within 30 arcmin of / = 
326.5°, b = 0°. 
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there is a 74 % chance that there will also be a 
12.2 GHz methanol maser, this 12.2 GHz detec- 
tion rate falls to 43 % for EGOs with associated 
6.7 GHz methanol masers but without OH masers. 
In comparison, co mparing the location s of EGOs 
with water masers iBreen et al.l (|2010bl ) has shown 
that 35 of the 89 wat er-methanol-OH sources t hat 
fell within the regions Cyganowski et al. ( 20081 ) in- 
spected for the presence of EGOs, equating to an 
association rate of 39 %. From these statistics it 
was concluded that EGOs persist into the stage 
of star formation that is evolved enough to have 
produced an OH maser, but not far past the stage 
where the developing UCHn region has caused the 
methanol maser emission to cease. Our results 
support this hypothesis along with the idea that 
12.2 GHz methanol masers are not present un- 
til well after the onset of the associated 6.7 GHz 
methanol maser emission. 

5. Conclusion 

We present the 12.2 GHz methanol maser 
properties of the 250 12.2 GHz methanol masers 
that we detect towards a complete sample of 580 
6.7 GHz methanol masers detected in the MMB 
survey. The detection rate of 43 % is lower than 
that of other large surveys of comparable sensitiv- 
ity but this discrepancy can be attributed to the 
biases introduced from the nature of the target 
sources in these searches. 

As found in previous searches, both the velocity 
ranges and the flux densities of the target 6.7 GHz 
sources surpass that of their 12.2 GHz companion 
in almost all cases. 80 % of the detected 12.2 GHz 
methanol maser peaks are coincident in velocity 
with the 6.7 GHz maser peak and the peaks of the 
remaining 12.2 GHz sources are usually associated 
with a close secondary feature. There are a cou- 
ple of instances where the peak of one or more of 
the 12.2 GHz methanol maser features surpasses 
that of the corresponding 6.7 GHz methanol maser 
peak. These unusual sources will be discussed in 
forthcoming catalog papers where we will explore 
in detail the specific data. 

We find that there is a non-uniform 12.2 GHz 
methanol maser detection rate throughout the 
Galactic plane. In particular we find that the 
12.2 GHz detection rate is highest in the longi- 
tude range 335 to 340°. Furthermore the longitude 



range 270 to 305° has a lower 12.2 GHz methanol 
maser detection rate than the 315 to 350° range. 

Some evidence of clustering of sources exhibit- 
ing 12.2 GHz masers (and likewise those with 
no 12.2 GHz emission) is apparent. We suggest 
that clustering may be a result of local bursts 
of star formation that have resulted in a num- 
ber of nearby sources being at similar evolutionary 
stages. 

Our analysis supports the evo lutionary sce- 
nario present in IBreen et al. ( 2010al) , whereby the 
12.2 GHz sources are associated with a somewhat 
later evolutionary stage than the 6.7 GHz sources 
devoid of this transition. Furthermore we find that 
both 6.7 GHz and 12.2 GHz methanol sources in- 
crease in luminosity as they evolve albeit more 
rapidly for the 6.7 GHz transition. In addition to 
this, evidence for an increase in velocity range as 
the sources evolve, is presented, and this also oc- 
curs at a greater rate for the 6.7 GHz sources. This 
implies that it is not only the luminosity, but also 
the volume of gas conducive to the different maser 
transitions, that increases as the sources evolve. 

Comparison with GLIMPSE sources has re- 
vealed a similar coincidence rate between the lo- 
cations of the masers and sources presented in the 
GLIMPSE point source catalog. It is evident that 
the GLIMPSE point sources that are associated 
with methanol masers have a large relative excess 
of 4.5 /im flux. We find little difference between 
the point source properties associated with the 
different maser categories (i.e. solitary 6.7 GHz 
sources, those with 12.2 GHz counterparts, and 
those associated with OH masers). There is a 
higher 12.2 GHz detection rate towards those 
6.7 GHz methanol masers that are coincident with 
EGOs. 

The Parkes radio telescope is part of the Aus- 
tralia Telescope National Facility which is funded 
by the Commonwealth of Australia for operation 
as a National Facility managed by CSIRO. This 
research has made use of NASA's Astrophysics 
Data System Abstract Service, the NASA/IPAC 
Infrared Science Archive (which is operated by 
the Jet Propulsion Laboratory, California Insti- 
tute of Technology, under contract with the Na- 
tional Aeronautics and Space Administration) and 
data products from the GLIMPSE survey, which 
is a legacy science programme of the Spitzer Space 
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